METHOD AND APPARATUS FOR REMOVING NOISE 
FROM FEATURE VECTORS 

BACKGROUND OF THE INVENTION 
The present invention relates to estimating 
the feature vectors corresponding to different 
sources that were combined to produce input feature 
vectors . 

A pattern recognition system, such as a 
speech recognition system, takes an input signal and 
attempts to decode the signal to find a pattern 
represented by the signal. For example, in a speech 
recognition system, a speech signal is received by 
the recognition system and is decoded to identify a 
string of words represented by the speech signal. 

To decode the incoming signal, most 
recognition systems utilize one or more models that 
describe the likelihood that a portion of the test 
signal represents a particular pattern. Typically, 
these models do not operate directly on the incoming 
signal, but instead operate on a feature vector 
representation of the incoming signal. In speech 
recognition, such feature vectors can be produced 
through techniques such as linear predictive coding 
(LPC) , LPC derived cepstrum, perceptive linear 
prediction (PLP), and mel-f requency cepstrum 
coefficients (MFCC) feature extraction. 

The incoming signal is often a combination 
of signals from different sources, each modified by a 
channel. For example, the incoming signal may be a 
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mixture of an original signal, which contains the 
pattern to be recognized, and one or more obscuring 
signals, such as additive noise and channel 
distortion. In speech recognition, the incoming 
5 signal may be a combination of the speech signal to 
be fed into a speech recognizer, additive noise, and 
channel distortion such as telephone channel 
distortion, or reverberations generated by the speech 
signal bouncing off walls in a room. Or, the incoming 
10 signal may be a combination of a speech signal with a 
Q channel signal (impulse response of the channel), 

JS where the channel signal is to be fed into a system 

H that recognizes channel types. Or, the incoming 

m signal may be a mixture of the speech signals from 

! Jf 15 two different speakers, each modified by a different 

v ■ channel, and each of which is to be fed into a speech 

|7| recognizer. 

W . Because noise and channel distortion make 

3 it more difficult to recognize a pattern in the 

20 incoming signal, it is often desirable to remove the 
noise and the channel distortion before performing 
pattern recognition. However, removing noise and 
channel distortion from the incoming signal itself is 
computationally difficult because of the large amount 
25 of data that has to be processed. To overcome this 
problem, some prior art techniques have tried to 
remove noise from the feature vector representation 
of the incoming signal instead of the incoming signal 
itself because the feature vector representation is 
30 more compact than the incoming signal. 
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However, past techniques for removing noise 
from feature vectors have relied on point models for 
the noise and the channel distortion. In other 
words, the noise reduction techniques have assumed 
5 that one single feature vector can represent the 
noise and another single feature vector can represent 
the channel distortion. The point models may be 
adapted to a sequence of input features, but they are 
held constant across the sequence. Because the noise 

10 and channel distortion vary across the sequence of 
input features, techniques that use this 
approximation do not accurately remove the noise or 
channel distortion. 

Some prior art techniques for removing 

15 noise from feature vectors attempt to identify the 
most likely combination of a noise feature vector, a 
channel distortion feature vector, and an original 
signal feature vector that would have produced the 
noisy feature vector. To make this determination, 

20 the prior art relies on an approximation of the 
relationship between noise, channel distortion, 
original signals, and incoming signals. 

However, prior art systems do not take the 
error present in the approximation into account when 

25 identifying possible combinations of noise, channel 
distortion, and original signals based on the 
incoming signal. In addition, the form of the 
approximation is typically set once and then used to 
identify the best combination. If the form of the 

30 approximation is not accurate, the resulting 
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identified combination of noise, channel distortion, 
and original signal will be inaccurate- However, the 
prior art does not provide a means for adjusting the 
form of the approximation to improve the resulting 
5 identified combination . 

SUMMARY OF THE INVENTION 
A method and computer-readable medium are 
provided for identifying clean signal feature vectors 
from noisy signal feature vectors. One aspect of the 

10 invention includes using an iterative approach to 
identify the clean signal feature vector. Another 
aspect of the invention includes using the variance 
of a set of noise feature vectors and/or channel 
distortion feature vectors when identifying the clean 

15 signal feature vectors. Further aspects of the 
invention use mixtures of distributions of noise 
feature vectors and/or channel distortion feature 
vectors when identifying the clean signal feature 
vectors. Additional aspects of the invention include 

20 identifying a variance for the noisy signal feature 
vector and using the variance when identifying the 
clean signal feature vector. 

BRIEF DESCRIPTION OF THE DRAWINGS 
25 FIG. 1 is a block diagram of a general 

computing environment in which the present invention 
may be practiced. 

FIG. 2 is a block diagram of a mobile 
device in which the present invention may be 
30 practiced. 



FIG. 3 is a block diagram of a speech 
recognition system in which one embodiment of the 
present invention is practiced. 

FIG. 4 is a flow diagram of the noise 
reduction technique of one embodiment of the present 
invention. 

FIG. 5 is a flow diagram of initialization 
step 450 of FIG. 4. 

FIG. 6 is a flow diagram of iteration step 

454 of FIG. 4. 

FIG. 7 is a graph showing a prior 
observation distribution, an observation distribution 
and posterior distributions during the iteration of 
FIG. 6. 

FIG. 8 is a graph showing prior 
distributions and posterior distributions for a 
mixture of components as well as a final mean for the 
combined posterior probability. 

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

FIG. 1 illustrates an example of a suitable 
computing system environment 100 on which the 
invention may be implemented. The computing system 
environment 100 is only one example of a suitable 
computing environment and is not intended to suggest 
any limitation as to the scope of use or 
functionality of the invention. Neither should the 
computing environment 100 be interpreted as having 
any dependency or requirement relating to any one or 
combination of components illustrated in the 
exemplary operating environment 100. 
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The invention is operational with numerous 
other general purpose or special purpose computing 
system environments or configurations. Examples of 
well known computing systems, environments, and/or 
5 configurations that may be suitable for use with the 
invention include, but are not limited to, personal 
computers, server computers, hand-held or laptop 
devices, multiprocessor systems, microprocessor-based 
systems, set top boxes, programmable consumer 
10 electronics, network PCs, minicomputers, mainframe 
computers, distributed computing environments that 
include any of the above systems or devices, and the 
like. 

The invention may be described in the 

15 general context of computer-executable instructions, 
such as program modules, being executed by a 
computer. Generally, program modules include 

routines, programs, objects, components, data 
structures, etc. that perform particular tasks or 

20 implement particular abstract data types. The 
invention may also be practiced in distributed 
computing environments where tasks are performed by 
remote processing devices that are linked through a 
communications network. In a distributed computing 

25 environment, program modules may be located in both 
local and remote computer storage media including 
memory storage devices. 

With reference to FIG. 1, an exemplary 
system for implementing the invention includes a 

30 general-purpose computing device in the form of a 



computer 110. Components of computer 110 may 

include, but are not limited to, a processing unit 
120, a system memory 130, and a system bus 121 that 
couples various system components including the 
system memory to the processing unit 120. The system 
bus 121 may be any of several types of bus structures 
including a memory bus or memory controller, a 
peripheral bus, and a local bus using any of a 
variety of bus architectures. By way of example, and 
not limitation, such architectures include Industry 
Standard Architecture (ISA) bus, Micro Channel 
Architecture (MCA) bus, Enhanced ISA (EISA) bus, 
Video Electronics Standards Association (VESA) local 
bus, and Peripheral Component Interconnect (PCI) bus 
also known as Mezzanine bus. 

Computer 110 typically includes a variety 
of computer readable media. Computer readable media 
can be any available media that can be accessed by 
computer 110 and includes both volatile and 
nonvolatile media, removable and non-removable media. 
By way of example, and not limitation, computer 
readable media may comprise computer storage media 
and communication media. Computer storage media 
includes both volatile and nonvolatile, removable and 
non-removable media implemented in any method or 
technology for storage of information such as 
computer readable instructions, data structures, 
program modules or other data. Computer storage 
media includes, but is not limited to, RAM, ROM, 
EE PROM, flash memory or other memory technology, CD- 
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ROM, digital versatile disks (DVD) or other optical 
disk storage, magnetic cassettes, magnetic tape, 
magnetic disk storage or other magnetic storage 
devices, or any other medium which can be used to 
5 store the desired information and which can be 
accessed by computer 100. 

Communication media typically embodies 
computer readable instructions, data structures, 
program modules or other data in a modulated data 

10 signal such as a carrier wave or other transport 
mechanism and includes any information delivery 
media. The term "modulated data signal" means a 
signal that has one or more of its characteristics 
set or changed in such a manner as to encode 

15 information in the signal. By way of example, and 
not limitation, communication media includes wired 
media such as a wired network or direct-wired 
connection, and wireless media such as acoustic, FR, 
infrared and other wireless media. Combinations of 

20 any of the above should also be included within the 
scope of computer readable media. 

The system memory 130 includes computer 
storage media in the form of volatile and/or 
nonvolatile memory such as read only memory (ROM) 131 

25 and random access memory (RAM) 132. A basic 
input/output system 133 (BIOS), containing the basic 
routines that help to transfer information between 
elements within computer 110, such as during start- 
up, is typically stored in ROM 131. RAM 132 

30 typically contains data and/or program modules that 
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are immediately accessible to and/or presently being 
operated on by processing unit 120. By way o 
example, and not limitation, FIG. 1 illustrates 
operating system 134, application programs 135, other 
5 program modules 136, and program data 137. 

The computer 110 may also include other 
removable /non-removable volatile /nonvolatile computer 
storage media. By way of example only, FIG. 1 
illustrates a hard disk drive 141 that reads from or 
10 writes to non-removable, nonvolatile magnetic media, 
3 a magnetic disk drive 151 that reads from or writes 

5 to a removable, nonvolatile magnetic disk 152, and an 

fi optical disk drive 155 that reads from or writes to a 

|| removable, nonvolatile optical disk 156 such as a CD 

g 15 ROM or other optical media. Other removable/non- 

; removable, volatile/nonvolatile computer storage 

7j media that can be used in the exemplary operating 

if environment include, but are not limited to, magnetic 

□ tape cassettes, flash memory cards, digital versatile 

20 disks, digital video tape, solid state RAM, solid 
state ROM, and the like. The hard disk drive 141 is 
typically connected to the system bus 121 through a 
non-removable memory interface such as interface 140, 
and magnetic disk drive 151 and optical disk drive 
25 155 are typically connected to the system bus 121 by 
a removable memory interface, such as interface 150. 

The drives and their associated computer 
storage media discussed above and illustrated in FIG. 
1, provide storage of computer readable instructions, 
30 data structures, program modules and other data for 
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the computer 110. In FIG . 1, for example, hard disk 
drive 141 is illustrated as storing operating system 
144, application programs 145, other program modules 
146, and program data 147. Note that these 

5 components can either be the same as or different 
from operating system 134 , application programs 135, 
other program modules 136, and program data 137. 
Operating system 144, application programs 145, other 
program modules 146, and program data 147 are given 

10 different numbers here to illustrate that, at a 
minimum, they are different copies. 

A user may enter commands and information 
into the computer 110 through input devices such as a 
keyboard 162, a microphone 163, and a pointing device 

15 161, such as a mouse, trackball or touch pad. Other 
input devices (not shown) may include a joystick, 
game pad, satellite dish, scanner, or the like. 
These and other input devices are often connected to 
the processing unit 120 through a user input 

20 interface 160 that is coupled to the system bus, but 
may be connected by other interface and bus 
structures, such as a parallel port, game port or a 
universal serial bus (USB) . A monitor 191 or other 
type of display device is also connected to the 

25 system bus 121 via an interface, such as a video 
interface 190. In addition to the monitor, computers 
may also include other peripheral output devices such 
as speakers 197 and printer 196, which may be 
connected through an output peripheral interface 190. 
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The computer 110 may operate in a networked 
environment using logical connections to one or more 
remote computers, such as a remote computer 180. The 
remote computer 180 may be a personal computer, a 
5 hand-held device, a server, a router, a network PC, a 
peer device or other common network node, and 
typically includes many or all of the elements 
described above relative to the computer 110. The 
logical connections depicted in FIG. 1 include a 
10 local area network (LAN) 171 and a wide area network 
(WAN) 173, but may also include other networks. Such 
networking environments are commonplace in offices, 
enterprise-wide computer networks, intranets and the 
Internet . 

15 When used in a LAN networking environment, 

the computer 110 is connected to the LAN 171 through 
a network interface or adapter 170. When used in a 
WAN networking environment, the computer 110 
typically includes a modem 172 or other means for 

20 establishing communications over the WAN 173, such as 
the Internet. The modem 172, which may be internal 
or external, may be connected to the system bus 121 
via the user input interface 160, or other 
appropriate mechanism. In a networked environment, 

25 program modules depicted relative to the computer 
110, or portions thereof, may be stored in the remote 
memory storage device. By way of example, and not 
limitation, FIG. 1 illustrates remote application 
programs 185 as residing on remote computer 18 0 . It 

30 will be appreciated that the network connections 



-12- 

shown are exemplary and other means of establishing a 
communications link between the computers may be 
used. 

FIG. 2 is a block diagram of a mobile 
5 device 200, which is an exemplary computing 
environment. Mobile device 200 includes a 
microprocessor 202, memory 204, input/output (I/O) 
components 206, and a communication interface 208 for 
communicating with remote computers or other mobile 
10 devices. In one embodiment, the afore-mentioned 
Q components are coupled for communication with one 

m another over a suitable bus 210. 

P Memory 204 is implemented as non-volatile 

111 electronic memory such as random access memory (RAM) 

s j2 15 with a battery back-up module (not shown) such that 

s -_ r information stored in memory 204 is not lost when the 

y general power to mobile device 200 is shut down. A 

Li: portion of memory 204 is preferably allocated as 

Q addressable memory for program execution, while 

20 another portion of memory 204 is preferably used for 
storage, such as to simulate storage on a disk drive. 

Memory 204 includes an operating system 
212, application programs 214 as well as an object 
store 216. During operation, operating system 212 is 
25 preferably executed by processor 202 from memory 204. 
Operating system 212, in one preferred embodiment, is 
a WINDOWS® CE brand operating system commercially 
available from Microsoft Corporation. Operating 
system 212 is preferably designed for mobile devices, 
30 and implements database features that can be utilized 
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by applications 214 through a set of exposed 
application programming interfaces and methods. The 
objects in object store 216 are maintained by 
applications 214 and operating system 212, at least 
5 partially in response to calls to the exposed 
application programming interfaces and methods. 

Communication interface 208 represents 
numerous devices and technologies that allow mobile 
device 200 to send and receive information. The 

10 devices include wired and wireless modems, satellite 
receivers and broadcast tuners to name a few. Mobile 
device 200 can also be directly connected to a 
computer to exchange data therewith. In such cases , 
communication interface 208 can be an infrared 

15 transceiver or a serial or parallel communication 
connection, all of which are capable of transmitting 
streaming information. 

Input/output components 206 include a 
variety of input devices such as a touch-sensitive 

20 screen, buttons, rollers, and a microphone as well as 
a variety of output devices including an audio 
generator, a vibrating device, and a display. The 
devices listed above are by way of example and need 
not all be present on mobile device 200. In 

25 addition, other input/output devices may be attached 
to or found with mobile device 200 within the scope 
of the present invention. 

FIG. 3 provides a block diagram of hardware 
components and program modules found in the general 

30 computing environments of FIGS. 1 and 2 that are 



particularly relevant to an embodiment of the present 
invention used for speech recognition. In FIG. 3, an 
input speech signal from a speaker 400 pass through a 
channel 401 and together with additive noise 402 is 
converted into an electrical signal by a microphone 
404 , which is connected to an analog-to-digital (A- 
to-D) converter 406. 

A-to-D converter 406 converts the analog 
signal from microphone 404 into a series of digital 
values. In several embodiments, A-to-D converter 406 
samples the analog signal at 16 kHz and 16 bits per 
sample, thereby creating 32 kilobytes of speech data 
per second. 

The output of A-to-D converter 406 is 
provided to feature extractor 400, which extracts a 
feature from the digital speech signal. Examples of 
feature extraction modules include modules for 
performing Linear Predictive Coding (LPC) , LPC 
derived cepstrum, Perceptive Linear Prediction (PLP) , 
Auditory model feature extraction, and Mel-Frequency 
Cepstrum Coefficients (MFCC) feature extraction. 
Note that the invention is not limited to these 
feature extraction modules and that other modules may 
be used within the context of the present invention. 

The feature extraction module receives the 
stream of digital values from A-to-D 406 and produces 
a stream of feature vectors that are each associated 
with a frame of the speech signal. In many 

embodiments, the centers of the frames are separated 
by 10 milliseconds. 



The stream of feature vectors provided by 
A-to-D converter 406 represents a noisy speech signal 
which is the combination of a clean speech signal, 
additive noise and channel distortion. These noisy 
feature vectors are provided to a noise reduction 
module 408 of the present invention, which generates 
a stream of "clean" feature vectors from the noisy 
feature vectors. 

The stream of "clean" feature vectors 
produced by noise reduction module 408 is provided to 
a decoder 412, which identifies a most likely 
sequence of words based on the stream of "clean" 
feature vectors, a lexicon 414, a language model 416, 
and an acoustic model 418. 

In some embodiments, acoustic model 418 is 
a Hidden Markov Model consisting of a set of hidden 
states. Each linguistic unit represented by the 
model consists of a subset of these states. For 
example, in one embodiment, each phoneme is 
constructed of three interconnected states. Each 
state has an associated set of probability 
distributions that in combination allow efficient 
computation of the likelihoods against any arbitrary 
sequence of input feature vectors for each sequence 
of linguistic units (such as words) . The model also 
includes probabilities for transitioning between two 
neighboring model states as well as allowed 
transitions between states for particular linguistic 
units. By selecting the states that provide the 
highest combination of matching probabilities and 
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transition probabilities for the input feature 
vectors, the model is able to assign linguistic units 
to the speech. For example, if a phoneme was 
constructed of states 0, 1 and 2 and if the first 
5 three frames of speech matched state 0, the next two 
matched state 1 and the next three matched state 2, 
the model would assign the phoneme to these eight 
frames of speech. 

Note that the size of the linguistic units 
10 can be different for different embodiments of the 
Q present invention. For example, the linguistic units 

jjj may be senones, phonemes, noise phones, diphones, 

^ triphones, or other possibilities. 

i y 

jjf In other embodiments, acoustic model 418 is 

r $ 15 a segment model that indicates how likely it is that 

a sequence of feature vectors would be produced by a 
hi segment of a particular duration. The segment model 

L!f differs from the frame-based model because it uses 

Q multiple feature vectors at the same time to make a 

^ 20 determination about the likelihood of a particular 

segment. Because of this, it provides a better model 
of large-scale transitions in the speech signal . In 
addition, the segment model looks at multiple 
durations for each segment and determines a separate 
25 probability for each duration. As such, it provides 
a more accurate model for segments that have longer 
durations. Several types of segment models may be 
used with the present invention including 
probabilistic-trajectory segmental Hidden Markov 
30 Models. 
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Language model 416 provides a set of 
likelihoods that a particular sequence of words will 
appear in the language of interest. In many 

embodiments , the language model is based on a text 
5 database such as the North American Business News 
(NAB) , which is described in greater detail in a 
publication entitled CSR-III Text Language Model, 
University of Penn., 1994. The language model may be 
a context-free grammar or a statistical N-gram model 
10 such as a trigram. In one embodiment, the language 
model is a compact trigram model that determines the 
probability of a sequence of words based on the 
combined probabilities of three-word segments of the 
sequence. 

15 Based on the acoustic model, the language 

model, and the lexicon, decoder 412 identifies a most 
likely sequence of words from all possible word 
sequences. The particular method used for decoding is 
not important to the present invention and any of 

20 several known methods for decoding may be used. 

The most probable sequence of hypothesis 
words is provided to a confidence measure module 420. 
Confidence measure module 420 identifies which words 
are most likely to have been improperly identified by 

25 the speech recognizer, based in part on a secondary 
frame-based acoustic model. Confidence measure 

module 420 then provides the sequence of hypothesis 
words to an output module 422 along with identifiers 
indicating which words may have been improperly 

30 identified. Those skilled in the art will recognize 



that confidence measure module 420 is not necessary 
for the practice of the present invention. 

Although the noise reduction technique of 
the present invention can be used in noise reduction 
module 408 , the invention is not limited to being 
used in speech recognition. Those skilled in the art 
will recognize that the invention may be used in any 
appropriate pattern recognition system. In addition, 
the invention may be used during the training of a 
pattern recognizer as well as during the detection of 
patterns. Those skilled in the art will also 
recognize that the method can be extended to multiple 
sources and multiple channels. Also, the invention 
may be used for purposes other than automatic 
recognition, such as denoising features for the 
purpose of reconstructing an enhanced signal. 

The noise reduction technique of the 
present invention identifies an optimal combination 
of obscuring signals and clean speech given an 
observed noisy speech vector. To do this, the 
present invention relies in part on the well-known 
Baysian rule: 

/ , \ p{y \n,x,c)p{n,x,c) 
p{n 9 x 9 c\y)= pr EQ. 1 

p\y) 

where p{n,x,c\y) is the posterior probability of a clean 
feature vector, x, a noise feature vector, n, and a 
channel distortion feature vector, c, given a noisy 
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feature vector, y, p{y\n,x,c) is the observation 

probability of a noisy feature vector given a clean 
feature vector, a noise feature vector, and a channel 
distortion feature vector, p(n 9 x 9 c) is the prior 
5 probability of the combination of a clean feature 
vector, a noise feature vector, and a channel 
distortion feature vector, and p(y) is the total 
probability for an observed noisy feature vector. 
Note that in this example, the noise feature vectors 

10 and the channel distortion feature vectors each 
provide separate forms of obscuring feature vectors. 

Once an approximation to the posterior 
p(n 9 x 9 c\y) has been found, the maximum a posteriori 
features may be chosen . Or, in another embodiment , 

15 the average values of the features may be chosen. 
Those skilled in the art will recognize that other 
statistics can be extracted from the posterior. In 
fact, a representation of the posterior probability 
distribution itself can be fed into a recognition 

20 system. 

Since p(y) is the same for all combinations 
of noise feature vectors, clean signal feature 
vectors, and channel distortion feature vectors, it 
can be ignored when searching for an approximation to 
25 the posterior distribution over noise feature 
vectors, clean signal feature vectors, and channel 
distortion feature vectors. 
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Under one embodiment of the present 
invention, the posterior probability (as well as the 
prior probability) is represented by a mixture of 
Gaussians, where each mixture component of the 
5 posterior probability is determined separately using 
a separate prior probability mixture component and a 
separate observation probability mixture component. 
Thus, mixture component i of the posterior 
probability is formed based on mixture component i of 

10 the prior probability and mixture component i of the 
observation probability. 

The process for identifying the combination 
of noise, channel distortion and original signal that 
provides the most likely posterior probability, 

15 p(n,x 9 c\y) is shown in FIG. 4. The process of FIG. 4 
begins at step 450 where the means and variances for 
the mixture components of the prior probability, 
observation probability and posterior probability are 
initialized. The process of step 450 is shown in 

20 greater detail in FIG. 5. 

In step 500 of FIG. 5, the means and 
variances for each mixture component of the prior 
probability are generated. To generate the means and 
variances, the process of the present invention first 

25 generates a mixture of Gaussians that describes the 
distribution of a set of training noise feature 
vectors, a second mixture of Gaussians that describes 
a distribution of a set of training channel 
distortion feature vectors, and a third mixture of 



Gaussians that describes a distribution of a set of 
training clean signal feature vectors. The mixture 
components can be formed by grouping feature vectors 
using a maximum likelihood training technique or by 
grouping feature vectors that represent a temporal 
section of a signal together. Those skilled in the 
art will recognize that other techniques for grouping 
the feature vectors into mixture components may be 
used and that the two techniques listed above are 
only provided as examples. 

After the feature vectors have been grouped 
into their respective mixture components, the mean 
and variance of the feature vectors within each 
component is determined. In an embodiment in which 
maximum likelihood training is used to group the 
feature vectors, the means and variances are provided 
as by-products of grouping the feature vectors into 
the mixture components. 

After the means and variances have been 
determined for the mixture components of the noise 
feature vectors, clean signal feature vectors, and 
channel feature vectors, these mixture components are 
combined to form a mixture of Gaussians that 
describes the total prior probability. Using one 
technique, the mixture of Gaussians for the total 
prior probability will be formed at the intersection 
of the mixture components of the noise feature 
vectors, clean signal feature vectors, and channel 
distortion feature vectors. 
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Once the means and variances for the 
mixture components of the prior probability have been 
determined, the process of FIG. 5 continues at step 
502 where initial means for the mixture components of 
5 the posterior probability are set. Under one 

embodiment of the invention, the initial means are 
set to be equal to the means of the prior 
probability T s mixture components. 

At step 504, the variances for the mixture 
10 components of the observation probability are 
determined. Under one embodiment, these variances 
are formed using a closed form expression of the 
form: 

VAR(y\x,n) = - EQ. 2 

15 where a is estimated from the training data. 

Under other embodiments, these variances 
are formed using a training clean signal, a training 
noise signal, and a set of training channel 
distortion vectors that represent the channel 

20 distortion that will be applied to the clean signal 
and noise signal. 

The training clean signal and the training 
noise signal are separately converted into sequences 
of feature vectors . These feature vectors, together 

25 with the channel distortion feature vectors are then 
applied to an equation that approximates the 
relationship between observed noisy vectors and clean 
signal vectors, noise vectors, and channel distortion 
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vectors. Under one embodiment, this equation is of 
the form: 

y *c+x+C (in (l + e( c "^-^ ]) )) EQ. 3 

where y is an observed noisy feature vector, c is a 
channel distortion feature vector, x is a clean 
signal feature vector, n is a noise feature vector, C 
is a transformation matrix, and C" 1 is the inverse of 
the transformation matrix. In equation 3: 



ln ( 1 + ,(c-b-^])) = 



In 



In 



A" 



In 



EQ. 4 



10 where n 3 , Cj, and Xj are the jth elements in the noise 
feature vector, channel feature vector, and clean 
signal feature vector, respectively, and C" 1 ^ is the 
i,j th entry of the inverse matrix C" 1 . 

In one embodiment of equation 3 above, the 

15 transform matrix C is an orthonormal matrix of 
discreet cosine transformation coefficients when the 
feature extraction technique produces cepstrum 
feature vectors. For embodiments that use a log 
spectrum feature extraction technique, C is the 

20 identity matrix. Those skilled in the art will 
recognize that other transform matrices, C, will be 
used in equation 3 depending on the particular 
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feature extraction technique used to form the feature 
vectors . 

In fact, C" 1 need not be square or of full 
rank, in which case C" 1 is a pseudoinverse matrix or 
5 another suitable matrix. 

Under one embodiment, the training clean 
signal feature vectors, training noise feature 
vectors, and channel distortion feature vectors used 
to determine the mixture components of the prior 
10 probability, are reused in equation 3 to produce 
C| calculated noisy feature vectors. Thus, each mixture 

component of the prior probability produces its own 
H set of calculated noisy feature vectors. 

m The training clean signal is also allowed 

*Z 15 to pass through a training channel before being 

si combined with the training noise signal- The 

jTj resulting analog signal is then converted into 

!'U feature vectors to produce a sequence of observed 

g noisy feature vectors. The observed noisy feature 

^ 20 vectors are aligned with their respective calculated 

noisy feature vectors so that the observed values can 
be compared to the calculated values. 

For each mixture component in the prior 
probability, the average difference between the 
25 calculated noisy feature vectors associated with that 
mixture component and the observed noisy feature 
vectors is determined. This average value is used as 
the variance for the corresponding mixture component 
of the observation probability. Thus, the calculated 
30 noisy feature vector produced from the third mixture 



component of the prior probability would be used to 
produce a variance for the third mixture component of 
the observation probability. At the end of step 504, 
a variance has been calculated for each mixture 
component of the observation probability. 

After the mixture components of the prior 
probability, observation probability, and posterior 
probability have been initialized, the process of 
FIG. 4 continues at step 452 where the first mixture 
component of the prior probability and the 
observation probability is selected. 

At step 454, an iteration is performed to 

find the mean for the posterior probability, p{n 9 x,c\y) 
of the selected mixture. The process for performing 
this iteration is shown in FIG. 6. 

In step 600 of FIG. 6, the prior 
probability, the observation probability, and the 
last determined mean of the posterior probability are 
used to identify a new mean for the posterior 
probability. In particular, the new mean for the 
posterior probability is calculated according to the 
variational inference principle and procedures as: 

EQ. 5 

where 77 is the newly calculated mean for the 
posterior probability of the current mixture, r] p is 
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the past mean for the posterior probability, X is 
the inverse of the covariance matrix for this mixture 
component of the prior probability, ju is the mean for 
this mixture component of the prior probability, ¥ is 
the variance of this mixture component of the 
observation probability, g(v p ) is the right-hand side 

of equation 3 evaluated at the last mean, g'^7 p ) is the 

matrix derivative of equation 3 calculated at the 
last mean, and y is the observed feature vector. 

In equation 5, ju , tj and 7j p are M-by-1 

matrices where M is three times the number of 
elements in each feature vector. In particular, ju , tj 

and 7j p are described by vectors having the form: 



-Elements For Clean Signal Feature Vector 
M 

— Elements For Noise Feature Vector 
3 

—Elements For Channel Distortion Feature Vector 
3 



EQ. 6 



Using this definition for ju , 77 and rj p , 
equation 3 above can be described as: 

M^H.- 1 — +l:M . - . - 



f 1 M 

+ H 1: T. 



+ Cln 



l+e 



V 



EQ. 7 
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where the designations in equation 7 indicate the 
spans of rows which form the feature vectors for 
those elements. 

In equation 5, the derivative g'iflp) ^ s a 

M 

matrix of order by-M where the element of row i, 

3 

column j is defined as: 



EQ. 8 



where the expression on the right side of equation 8 
is a partial derivative of the equation that 
10 describes the ith element of giflp) relative to the jth 

element of the r/ p matrix. Thus, if the jth element 

of the rf p matrix is the fifth element of the noise 

feature vector, n 5 , the partial derivative will be 
taken relative to n 5 . 
15 Note that when the transform matrix, C, of 

equation 7 is equal to the identity matrix, the ith 
element of giflp) is defined as: 



20 



g 



( \ 



=c, +x 



, +ln(l 



EQ. 9 



so that the partial derivative only has nonzero 
values for [t] p ]j equal to ru, Ci, or Xi. 
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After equation 7 has been used to determine 
a new mean for the posterior probability at step 600 , 
the process of FIG. 6 continues at step 602 where a 
stopping criterion is tested. For example, the mean 
5 may be examined to determine whether it has 
converged. Under most embodiments, this is 

determined by comparing the new mean to the last mean 
determined for the posterior probability. If the 
difference between these two means is larger than 

10 some threshold, the process returns to step 600 using 
the new mean as the last mean and determining a 
revised new mean. In another embodiment, a fixed 
number of iterations is performed. Those skilled in 
the art will recognize that other techniques for 

15 determining when to stop the iterations may be used 
with the present invention. 

Steps 600 and 602 are repeated until the 
stopping criterion is satisfied, at which point the 
iteration of FIG. 6 ends at step 604. When the 

20 process of FIG. 6 reaches step 604, step 454 of FIG. 
4 is complete. The process of FIG. 4 then continues 
at step 456 where the variance for the posterior 
probability of the current mixture component is 
determined. Under one embodiment of the present 

25 invention, the variance for the posterior probability 
is found in the premultiplier in the second term of 
the right hand expression of equation 5. The 
variance is determined by evaluating this factor at 
the selected mean for the posterior probability. 
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The effects of the iterations of FIG. 6 are 
shown in FIG. 7 for a single feature of the feature 
vectors. In FIG. 7, feature i of a clean signal 
feature vector is shown along horizontal axis 700 
5 while feature i of a noise feature vector is shown 
along vertical axis 702. A distribution 704 with a 
mean 7 06 is shown in this space for a mixture 
component of the prior probability. FIG. 7 also 
includes a distribution 705 for the observation 

10 probability. 

Before the iterations of FIG. 6 begin, the 
mean of the posterior probability is set equal to 
mean 706. After a first iteration through step 600, 
the mean of the posterior probability has shifted to 

15 location 708. After a second iteration, the mean has 
shifted to location 710. On the final iteration, the 
mean moves to location 712 and the iteration ends. 
The variance of the posterior probability is then 
determined at step 456 providing a distribution 714 

20 for the posterior probability. 

Note that the mean of the posterior 
probability settles at a location between the prior 
probability distribution and the observation 
probability distribution 705. Thus, the posterior 

25 probability distribution is a balance between the 
prior probability distribution 704 and the 
observation probability 705. 

After the mean and variance for the first 
mixture component of the posterior probability has 

30 been determined, the process of FIG. 4 continues by 
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determining whether there are more mixture components 
at step 458. If there are more mixture components, 
the next mixture component is selected at step 4 60 
and steps 454 and 456 are repeated for the new 
mixture component . 

Once all of the mixture components have had 
their mean and variance determined for the posterior 
probability, the process of FIG. 4 continues at step 
4 62 where the mixture components are combined to 
identify a most likely noise feature vector, channel 
distortion feature vector, and clean signal feature 
vector given the observed noisy signal feature 
vector. Under one embodiment, the noise feature 
vector, channel feature vector, and clean signal 
feature vector are calculated as: 



"post = ^Ps!h 



(M , 2M 

— +1: 

I 3 3 . 



EQ. 10 



EQ. 11 



20 



post 



5=1 



EQ. 12 



where S is the number of mixture components, p s is 



the weight for mixture component s, ij. 



M , 2M 

— +1: 

3 3 



is 



25 the noise feature vector for the mean of the 
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'2M 



is the channel 



posterior probability, tj^ — - + l:M 

distortion feature vector for the mean of the 

M 



posterior probability, t/J l: 



is the feature vector 



for the clean signal of the mean of the posterior 
5 probability, and n post , c pos t, and x pos t are the average 
values of the noise feature vector, channel 
distortion feature vector, and clean signal feature 
vector, respectively, given the observed noisy 
feature vector. 

10 The weight for each mixture component, p s 

is calculated as: 



7te Gs 



/=1 



15 where the dominator of equation 13 normalizes the 

weights by dividing each weight by the sum of all 

other weights for the mixture components. In 

equation 13, 7t s is a weight associated with the 
mixture components of the prior probability and is 

20 determined as: 

n. =71* ■< EQ. 14 



where n* , n" , and < are mixture component weights for 
the prior clean signal, prior noise, and prior 
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channel distortion, respectively. These weights are 
determined as part of the calculation of the mean and 
variance for the prior probability. 



5 affects the weighting of a mixture component based on 
the shape of the prior probability and posterior 
probability, as well as the similarity between the 
selected mean for the posterior probability and the 
observed noisy vector and the similarity between the 
10 selected mean and the mean of the prior probability. 

Under one embodiment, the expression for G s is: 



In equation 13, 



is a function that 



2 2 1 




15 





20 



1 



if diagonal element of (g (rj^f *¥~ ] g (^)s; 1 



sum o\ 



2 



EQ. 15 
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where ln|2^Z,| involves taking the natural log of the 
determinate of 2n times the covariance of the prior 
probability, ln|2^<D s | involves taking the natural log 

of the determinant of 2k times the covariance matrix 
5 of the posterior probability, which is the 
premultiplier of the second term of the right hand 
side of equation 5. 

Those skilled in the art will recognize 
that there are other ways of using the mixture 

10 approximation to the posterior to obtain statistics. 
For example, the means of the mixture component with 
largest p can be selected. Or, the entire mixture 
distribution can be used as input to a recognizer. 

An example of the determination of the 

15 posterior noise feature vector, channel distortion 
feature vector, and clean signal feature vector are 
shown in FIG. 8. In FIG. 8, as in FIG. 7 above, 
feature i of the clean signal is shown along 
horizontal axis 800 and feature i of the noise is 

20 shown along vertical axis 802. Note that for 
simplicity, feature i for the channel is not shown, 
but would provide a third dimension if placed in FIG. 
8. 

In FIG. 8, there are six mixture components 
25 804, 806, 808, 810, 812 and 814 for the prior 
probability. The prior mixture components are 

associated with six mixture components for the 
posterior probability indicated as distributions 816, 
818, 820, 822, 824 and 826, respectively. These 
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posterior mixture probabilities are combined to 
identify a single vector 828 that describes the most 
likely clean signal feature vector, noise feature 
vector, and channel distortion feature vector given 
5 the observed noisy feature vector. Note that in FIG. 
8, only one feature is shown, however the discussion 
of FIG. 8 should be interpreted as extending to all 
of the features of the feature vectors . Thus, in 
practice, FIG. 8 is an M dimensional space and vector 

10 828 is an M dimensional vector. 

After the noise feature vector, channel 
distortion feature vector, and clean signal feature 
vectors have been determined from the posterior 
mixture components, the process of FIG. 4 continues 

15 at step 4 64 by determining if there are anymore noisy 
vectors that need to be cleaned. If there are, steps 
452, 454, 456, 458, 460 and 462 are repeated to 
generate a clean signal vector for the noisy vector. 
When all of the noisy vectors have been processed at 

20 step 464, the noise reduction technique of FIG. 4 
ends at step 466. 

The present invention differs from the 
prior art in a number of ways. First, as discussed 
above, the present invention utilizes the variance of 

25 the noise and channel features. Also, the present 
invention utilizes a mixture of Gaussians to 
represent the noise component of the prior 
probability and to represent the channel distortion 
component of the prior probability. By using 

30 mixtures of Gaussians to model the noise and channel 



distortion, it is thought that the present invention 
will remove noise and channel distortion more 
accurately than if the noise and channel were modeled 
as single points as was done in the prior art. In 
practice, it is highly likely that the use of mixture 
of Gaussians for noise and for channels allows the 
algorithm to deal effectively with time-varying noise 
and channels, because an instantaneous noise and 
channel value can be accurately represented by one of 
many Gaussian components in the mixture distribution. 

Similarly, the present invention, as shown 
in equation 5, also takes the variance of the 
observation probability into consideration when 
identifying the mean for the posterior probability. 
In prior art noise reduction techniques, this 
variance was assumed to be zero. By taking this 
variance into account, the present invention is able 
to more accurately select the mean for the posterior 
probability because it takes into account the error 
present in the approximation of equation 3. 

Lastly, the iterative technique of the 
present invention is not shown in the prior art. 
Thus, prior art noise reduction systems do not 
iteratively modify the estimate of the clean signal 
vector. Instead, the prior art makes a single 
selection for the clean feature vector and does not 
try to improve upon that selection once it has been 
made . 

Although the invention has been described 
above with reference to two signals (a clean signal 



and a noise signal), and one channel, the invention 
is not limited to this combination. In particular, 
additional signals from additional sources may be 
present and the signals may pass through more than 
one filter or channel. Those skilled in the art will 
recognize that the equations described above may be 
extended to cover any number of signals and any 
number of channels. 

Although the present invention has been 
described with reference to particular embodiments, 
workers skilled in the art will recognize that 
changes may be made in form and detail without 
departing from the spirit and scope of the invention. 



